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Abstract. In recent years inelastic spin-flip spectroscopy using a iow- 
temperature scanning tunneiing microscope has been a very successful tool for 
studying not only individual spins but also complex coupled systems. When these 
systems interact with the electrons of the supporting substrate correlated many- 
particle states can emerge, making them ideal prototypical quantum systems. 
The spin systems, which can be constructed by arranging individual atoms 
on appropriate surfaces or embedded in synthesized molecular structures, can 
reveal very rich spectral features. Up to now the spectral complexity has 
only been partly described. This manuscript shows that perturbation theory 
enables one to describe the tunneling transport, reproducing the differential 
conductance with surprisingly high accuracy. Well established scattering models, 
which include Kondo-like spin-spin and potential interactions, are expanded to 
enable calculation of arbitrary complex spin systems in reasonable time scale and 
the extraction of important physical properties. The emergence of correlations 
between spins and, in particular, between the localized spins and the supporting 
bath electrons are discussed and related to experimentally tunable parameters. 
These results might stimulate new experiments by providing experimentalists with 
an easily applicable modeling tool. 
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1. Introduction: 

With low-temperature scanning tunneling microscopes (STM) scientist have developed 
a tool that has the ability to detect and manipulate individual magnetic spin systems 
on the atomic and molecular level by an externally controlled probe with sub-nm 
precision. These instruments have opened a new field of research envisioning not 
only a deeper understanding of the origin of molecular magnetism by studying the 
interactions between nanoscale spins but also of the many-particle effects between the 
localized spins and the itinerant electrons of the supporting substrate. The progress in 
this field is best demonstrated by the recently achieved ability to build stable magnetic 
bits by cleverly arranging only a handful of Fe atoms on either a thin insulating [T] or 
on a nonmagnetic Cu(lll) surface [2]. While in these experiments the transition from 
quantum-mechanical to classical magnetic behavior is explored [3] , other applications 
might be envisioned ranging from spin-based logic circuits [4] to entangled systems in 
which the quantum mechanical nature is crucial for computational purposes [5]. 

At the basis of all these experiments lies the spectroscopic capability of the 
STM. About ten years ago A. Heinrich and his co-workers showed in a hallmark 
experiment that it is possible to detect inelastic spin-flip excitations on Mn atoms 
adsorbed on patches of AI 2 O 3 on a NiAl surface [5]. Since then many experiments 
have focused on transition metal atoms on a thin layer of CU 2 N on Cu(lOO) 
[IlIHlIIllinillllllllllllllllllllllllllZlIIHlIIlllO]- In these experiments, patches of 
CU 2 N are formed by sputtering a clean Cu(lOO) surface with N+ ions and subsequent 
annealing. This leads to a monolayer of CU 2 N on the surface on which the desired metal 
atoms are deposited (Figure [T^). Experimentally, every 3d transition metal atom 
adsorbed on this surface reveals its fingerprint when a spectrum is taken by placing 
the tunneling tip of an STM over the atom (Figure [lb). The spectra are measured by 
varying the bias voltage applied between tip and sample and recording the differential 
conductance dl/dV. All spectra show characteristic features that can be classified 
into two groups: (1) Step-like increases in the differential conductance, positioned 
symmetrically around zero bias, and ( 2 ) peaks of the differential conductance at zero 
bias. 

The number of distinguishable conductance steps varies for different spin systems 
and the precise form of the steps often shows some asymmetry with respect to the 
applied bias direction and can exhibit some overshoot of the conductance at the 
step-energy. As we will see in detail, the steps are due to the opening of additional 
conductance channels precisely governed by the magnetic properties of the spin system. 
These excitations can be present even at zero magnetic field due to the magneto¬ 
crystalline anisotropy. The anisotropy is caused by the reduction of the geometric 
symmetry at the surface and by spin-orbit coupling [23l [24l [25], which have the 
effect of lifting the inherent degeneracy of the spin states. For single atoms the 
maximum anisotropy is limited to a few tens of meV [26], whereby the magnetic 
excitation energies usually range from less than one meV up to a few meV requiring 
experimental setups operating at temperatures T < 4 K. Here, it is worth to note 
that one has to be aware that not all energetically low lying step-like increases in 
the differential conductance must originate from magnetic excitations. The tunneling 
electrons can also excite low-energy mechanical vibrations which can produce similar 
strong spectroscopic features llzlEHlElllSol HaisT] and which can interact with 
the spin excitations [321 133 ] • Thus, to clearly distinguish magnetic excitations their 
behavior in external applied magnetic fields is often crucial. 
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Figure 1. Some examples of inelastic spin-flip spectroscopy on single atoms and 
molecules, (a) Constant current topography of the Cu 2 N surface on Cu(lOO) 
with co-adsorbed 3d metal atoms, (size approx.: 15 x 40 nm^, V = 10 mV, 
7 = 1 nA). (b) Differential conductance dljdV spectra measured on top of 
Mn, Fe, Ti, and Co adatoms (from top to bottom) at zero magnetic field and 
T = 0.6 K. Each atom reveals a characteristic fingerprint in the dljdV signal, 
(c) Constant current topography of a Mni 2 Aci 6 molecular magnet adsorbed on 
h-BN on Rh(lll) (V = -1 V, 7 = 45 pA). (d) Typical dl/dV and d?I/dV‘^ 
spectra at T = 1.5 K and zero magnetic field on Mni 2 Aci 6 . (e) Constant current 
topography of an organic radical molecule (C 28 H 25 O 2 N 4 ) adsorbed on Au(lll) 
{V = 100 mV, 7 = 33 pA, Contour line distances 50 pm).(f) dljdV spectra at 
different temperatures (T = 1.5, 3, 4.4, 6.7 K, bottom to top, upper panel) and 
magnetic fields {3 = 14, 10, 6 , 2 T, bottom to top, lower panel) on the radical. 
Spectra are vertically offset for clarity. Data adapted from [O] [ 8 ] 1211122| . 


The peaks at zero bias are due to the Kondo effect in which itinerant electrons 
from the substrate coherently scatter with the localized magnetic moment of the 
adatom [22 [351 [Ml. This effect has first been detected by scanning tunneling 
spectroscopy on single metal atoms adsorbed on noble metal substrates [MIIMIIMI 
EDI HU da dal mi sa [46]. In these early experiments the metal atoms are relatively 
strongly coupled to the substrate leading to a characteristic Kondo temperature in the 
range of Tk ~ 30 — 300 K, as determined by the full-width at half-maximum of the 
peak, and in most cases to a strongly asymmetric Fano lineshape due to interference 
effects with a potential scattering channel 071 Ei ED]. A decoupling layer, such 
as CU 2 N, significantly reduces the Kondo temperature suggesting the possibility of 
influencing this many-body state with experimentally accessible magnetic fields and 
temperatures [iEDllISlEI]. This enables the study of the interplay between Kondo 
screening, the magnetic anisotropy, and nearby spins uni HD]. 

Apart from CU 2 N, other decoupling substrates have been exploited to study spin 
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excitations. For example, low-lying excitations associated with a spin S' = 1 of the 
magnetic atom have been found for single Fe atoms embedded into the semiconducting 
InSb(llO) surface [HU [53], which exhibits a two-dimensional electron gas confined to 
the surface. Other experiments used two-dimensional materials, such as graphene 
and hexagonal boron nitride (/i-BN), as substrates for metal adatoms or magnetic 
molecules. For example, the Kondo state was observed for individual Co atoms 
adsorbed on a graphene sheet on top of a Ru(OOOl) surface with different Kondo 
temperatures and effective ^-factors depending on the adsorption site with respect 
to the underlying metal [SO]. In a different experiment Co and CoHa, (a; = 1 — 3) 
complexes on graphene on top of Pt(lll) have been investigated revealing a S' = I 
for Co and C 0 H 3 and a high magnetic anisotropy of ~ 8.1 meV for the Co adsorbate 
[51] . On the highly corrugated and insulating /i-BN adlayer on top of a Rh(lll) 
substrate, CoH and C 0 H 2 complexes showed both, spin-flip excitations and the Kondo 
effect, pointing to two different spin states of S = I and S = 1/2, respectively [5T] . 
Interestingly, the CoH complex revealed a dependence of the effective anisotropy on 
the coupling to the underlying substrate, similarly as observed for Co atoms on large 
patches of CU 2 N [ 16 ] . 

Spin excitations can also be observed for spin systems adsorbed on bare metal 
substrates, even though lifetime, anisotropy energies, and intensities are in general 
reduced due to the strong coupling with the substrate. For example, spin-excitations 
of Co and Fe on Pt(lll) have been measured [55] [56], where the measurements of 
reference [55] were performed at T « 6 K and showed an approximately ten-times 
higher apparent magnetic anisotropy as the latter measurement (reference [56] 1 that 
was performed at T = 0.3 K. This discrepancy stems from the strong temperature 
dependent broadening of the spectrum, which leads to this gross overestimation 
and illustrates the need for low-temperature measurements m- Additionally, spin 
excitations have been detected for Fe atoms adsorbed on Cu(lll) [58] and on Ag(III) 
[59] . Apart from the 3d transition metal atoms, 4f lanthanide atoms also show low- 
energy magnetic signals. For example, the Kondo effect was measured on small Ce 
clusters on CU 2 N m and spin-flip signals were observed for Gd adsorbed on Pt(III) 
and Cu(lll) [5D] and Ho on Pt(lll) [5T], even though the intensities of the spin-flip 
excitations are unclear but presumably very small due to the strong localization of 
the 4f wavefunctions close to the nuclei and their small spacial extension [62] . 

In addition to these, metal-organic complexes, such as M-phthalocyanine with 
M = Mn, Fe, Co, Ni, and Cu, have been studied on different surfaces [63l[64ll651l66ll67l 
[M]|M]IZD]- Many of them showed Kondo screening with a lower Kondo temperature 
enabling the tuning by dehydrogenation [63] EQ] or the splitting of the peak by an 
externally applied magnetic field [641168] or by coupling to a ferromagnetic substrate 
m- Interestingly, FePc adsorbed on Au(lll) showed a clear Kondo signature [69] . 
while the same molecule adsorbed on top of a CuO decoupling layer showed a double¬ 
step in the differential conductance pointing to S' = I [55] . 

The electronic decoupling can also be created by using a substrate in the 
superconducting state in which the creation of Cooper pairs leads to a strongly reduced 
density of unpaired electrons around the Fermi energy. This results in a significant 
increase of the spin lifetime of the metal-organic complex m- Furthermore, such a 
substrate allows one to study the competition between superconducting phenomena 
and Kondo screening as observed on MnPc adsorbed on Pb(III) [571175] . 

Apart from metal atoms and metal-organic complexes spin-flip excitations were 
detected in complex molecules in which several spin centers are coupled forming 
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a giant spin, as it has been shown at the example of the prototypical molecular 
magnet manganese-12-acetate-16 which has a ground state with a total spin of 
^ = 10 [5T] (Figure [TJ; and d). Here, the /i-BN decoupling layer is crucial 
as the magnetism is strongly quenched when the molecule is adsorbed upon a 
Au(lll) surface. Furthermore, fully organic molecules including the charge-transfer 
complexes TTF-TCNQ (tetrathiafulvalene-tetracyanoquinodimethane) [32] and an 
organic radical molecule [22], have shown the Kondo effect, whereby for the latter 
molecule the temperature and magnetic field dependence was fully understood and 
simulated in a perturbative approach (figure |T^ and f). 

To describe the inelastic spin-flip excitation spectra, scattering models have been 
developed that treat the interaction of the tunneling electron with the localized 
spin system effectively as a one-electron second-order perturbation [711131711177]. 
Additionally, similar models allowed rationalization of the change in the spectra when 
a spin-polarized tip is employed [75] and the dynamics at increased coupling between 
tip and sample [711|5ni[5I]. To address the experimentally observed bias asymmetry, 
co-tunneling models in the strong Coulomb-blockade regime have been proposed |82j 
and many-electron effects using the non-crossing approximation have been included 
[83]. Furthermore, third-order scattering models similar to the ones used here were 
employed to well describe the observed bias overshooting at certain inelastic excitation 
steps [511[5S]. Nevertheless, these models were restricted to Kondo-like interactions 
and did not include potential scattering. 

In this paper, I plan to review the straightforward model of the exchange- 
interaction between an isolated spin system and the tunneling electrons and use a 
perturbative tunneling approach to simulate experimental data with unprecedented 
accuracy. The basic idea of the model goes back 50 years to the hallmark discovery 
of Juan Kondo that higher order scattering of bulk electrons on a magnetic impurity 
leads to logarithmic divergences [Ml [35]. We will see that such a model enables 
the determination of physical properties like the magnetocrystalline anisotropy, 
the coupling strength to the substrate, the state lifetimes, and the magnetization 
directly by comparing the differential conductance measurements with simulations. 
Furthermore, it allows one to grasp some of the correlations and entanglements 
that are formed by the many-particle interactions due to the large electron bath of 
the substrate. The model, even though inherently limited due to its perturbative 
approach, has the advantage that we can develop it straightforwardly using simple 
matrix algebra. Due to this simplicity it is computationally very fast and can thus 
help us to understand future experimental observations. In these calculations the 
magnetic anisotropies, gyromagnetic factors, and coupling strengths to the substrate 
enter as experimentally determined parameters. Note however, that in particular 
for transition metal atoms adsorbed on metallic substrates the magnetic anisotropy 
have been successfully determined from first principle by, for example, time-dependent 
density functional theory [53 |57J [SSj [55] or density functional theory which includes 
spin-orbit coupling [59| [90] . 

The basic idea of the model is to treat the tunneling as a scattering event of an 
incoming electron with the localized spin system (Figure [2^). Postponing the exact 
meaning for later, it is amazing that only two parameters govern these interactions: 
The Kondo scattering parameter J and the potential scattering parameter U. These 
two parameters are crucial for the understanding of the low-energy excitations. One 
is tempted to ask “What is the origin of these?” To find an answer it is helpful to 
step back from the scattering picture and look at the single-impurity Anderson model. 
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Figure 2. Comparison of the Kondo scattering and the single-impurity Anderson 
model, (a) Schematic view of the Kondo model where electrons scatter via spin- 
flip or the potential interaction with a localized spin system. (b) Schematic view 
of the single orbital Anderson impurity model. An occupied 3d orbital located 
at an energy below the Fermi energy hybridize with the substrate electrons 
leading to a state broadening A. The occupation of this orbital with a second 
electron is prohibited by the Coulomb repulsion energy Ud- 


which describes one half-filled atomic orbital (for example from a 3d transition metal 
atom) interacting with the continuous states of a host metal (Figure [2 ]d) [HI]- In this 
picture all energies are in the eV range and it is not per se clear how this relates to 
the scattering model with low-energy excitations in the meV range. These relations 
were found by Schrieffer and Wolf, who were able to relate the high energies of the 
Anderson model with the scattering parameters |92) : 



1 


Ed + Ud 


< 0 , 



Ed + Ud 


( 1 ) 


Here it is crucial for the understanding of the forthcoming model that the Kondo 
spin-spin exchange scattering J is for 5" = 1/2 systems always negative, which means 
that the localized spin and the substrate electrons couple antiferromagnetically, while 
the potential scattering U can have both signs, either positive or negative. However, 
this result can be different for higher spins. 

The manuscript is organized as follows: In section [5] the basic model in the zero- 
current approximation is introduced and discussed for the example spin S = 1/2. 
In section [3] high-spin systems are discussed. Experimentally observed examples are 
given and compared to the model calculations. Section |4] discusses the limitations 
of the model, in particular, its inability to cover all correlation effects that occur 
when a system enters the strong Kondo regime using experimental data obtained on 
the C 0 /CU 2 N system as an example. In section [S] the initial model is extended by 
including rate equations and tested against experimental data. Here, we will observe 
that the tunneling current too can lead to the appearance of a non-equilibrium Kondo 
effect. Finally, section [6] summarizes the manuscript and outlines possible routes for 
future extensions. 


2. The model 

To describe the experimental observations we use a simplified model in which the 
Hamiltonian of the total quantum mechanical system is divided into the ones of the 
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subsystems of the two electron reservoirs in tip and sample, the one of the localized 
spin system, and an interaction Hamiltonian that enables the exchange of charge 
carriers between the reservoirs: 


H = Ht + Hs +Ha+ H' 


( 2 ) 


The Hamiltonian of the tip Ht and sample Hg, respectively, can be described in the 
framework of second quantization using creation and annihilation operators and a: 

Ht = ( 2 ) 

k,(T 

Hs = 'y ^ esfeo-dg^^asfco-, (4) 

k^cr 


with etfccr and Cska- as the energy of the electrons with momentum k and spin a in the 
tip and sample, respectively. 

These two many-particle systems are the source and sink for the tunneling 
electrons in the STM experiment. Instead of using creation and annihilation operators 
in the momentum space fc, we assume in the small energy range of interest, i. e. to 
some tens of meV around the Fermi energy, for tip and sample a continuous and 
energetically flat density of states p{e) = {a\^ak'a^ = p, with (•) as 

the time averaged expectation value. 

In general, the electronic states |p) in these electrodes might be spin-polarized in 
an arbitrary direction which we account for by the corresponding spin density matrices 
g = |p)(p| = \{I -I- fi ■ a). Here, n describes the direction and —1 < |n| < 1 the 
amplitude of the polarization in the chosen coordinate system, & = (ax, are the 

standard Pauli matrices, and I is the (2 x 2) identity matrix. With this convention the 
spin polarization is identical to the relative imbalance between majority and minority 
spin densities 77 = |n| = 


PT~Pr 

Pt+Pi 


where f and j, account for the two different spin 
directions along the chosen quantization axis. This description via density matrices 
allows arbitrary polarization directions in tip and sample which obey the quantum 
statistics (see [Appendix A] ). 

The impurity spin system may either contain only a single spin or a finite number 
of coupled spins. We describe this system by a model Hamiltonian Ha that includes 
Zeeman and magnetic anisotrcmy energy and - in the case of more than one spin - the 
Heisenberg coupling terms Jijjj and the non-collinear Dzyaloshinskii-Moriya coupling 

Hji 


'ly between individual spins: 


Ha = 2_^giPBB ■ 


E 


(^D,{sir+E,[{sir - {Sir]) 




(5) 


In this equation is the Bohr magneton and Di, and Ei are the gyromagnetic 
factor, the axial and the transversal m^netic anisotropy for the z-th spin in the 
reference coordinate frame, respectively^ The externally applied magnetic field is 


f We write as vector to enable also anisotropic Heisenberg or Ising-like couplings. 

I Note, that instead of using the rather phenomenological D and E values the model can be easily 
adapted to more physical operators that connect to the spin-orbit couplings m or to the symmetries 
of the system m- 









B and the total spin operators S* = {Sl.,Sy,Sl) are built from operators of the 
form Sl.y^^ = Ii ® ® ® ••• ® In, which only act on the i-th spin and 

where / denotes the identity matrix and 0 the Kronecker matrix product. The 
spin operators SI. y ^ can be easily constructed remembering that Szl^pm) = Ii-'mz\'<pm) 
with niz = S, S — 1, ■ ■ ■ ,—S as the magnetic quantum number along our chosen z- 
axis and S±\^rn) = ^\/{S =F m){m ± S + l)|'0m±i)- These then enable calculation of 
Sx = 5 (<5+ + S-) and Sy = —Ii{S+ — S-), respectively. For simplicity we will set 
h = 1 from now on. 

Diagonalizing the Hamiltonian of the localized spin system (equation [5]) leads to 
discrete energy eigenvalues e„ and eigenstates \ip)n- For a single spin there are 2S' + 1 
eigenstates, while for complex spin structures that consist of several spins the number 
of eigenstates increases quickly as + 1). Because we neglect direct interactions 

between the electron baths in tip, sample and the localized spin system, we describe 
the total state as a product of the continuous electron states \'.p) and the discrete spin 
states IV'), i-c. Note, however, that the coupling of the spin- 

system with the substrate will lead to an entanglement between sample electrons and 
the spin system (see section U]) and to a renormalization of the parameters in equation 
[5][9lliaiMl[5l]. Here we assume that the renormalization is already included in the 
anisotropies and gyromagnetic factors and omit for the moment the entanglement. 

The interaction Hamiltonian H' of equation allows for tunneling of electrons 
from tip to sample or vice versa only via Kondo-like spin-flip or potential scattering 
processes with the impurity (Figure ISt-): 

H' = + ^s->t + hs-s-s, 

aL atA , (6) 

4x' asA • s* + , (7) 

i,X',X ^ 

with Tl^ as the coupling constants between tip and the i-th adsorbate spin, and 
Ui = UijJi as the unitless ratio between Kondo and potential scattering. While 
we assume that the impurity is much more strongly coupled to the sample than to the 
tip, spin-spin scattering between the impurity and substrate electrons are additionally 
considered: 

= “>5A' a-■S\ (8) 

i,X',X 

With the above assumptions the model is largely identical to the ones studied 
by Appelbaum and Anderson already in the 70’s for mesoscopic tunnel junctions 

[ 131 ^ 117 ]. 

2.1. Current and conductance in second order 

We will now briefly review the calculation of the tunneling current using Fermi’s golden 
rule. In STM experiments the bias eV applied between tip and sample shifts the Fermi 
level CF of the tip with respect to the one of the sample. At positive electrons from 
occupied tip states can cross the tunnel-barrier and interact with the localized spin 
system under exchange of angular momentum and energy (Figure[3})). To obey energy 
conservation, the energy difference e^/ = e/ — between initial and final state of the 
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tipp sample p tipp sample p 


Figure 3. Model of the tunneling process between tip and sample, (a) Electrons 
originating from the tip interact with the localized spin system either via spin-spin 
(blue arrow) or potential scattering (red arrow). Additionally, we allow for spin- 
spin exchange scattering between sample electrons and the localized spin (yellow 
arrow), (b) Scheme of the first order Born approximation: In the tunneling 
process of an electron from the initial state i to the final state / it exchanges 
angular momentum and energy with the localized spin system, (c) In second 
order Born approximation the system additionally occupies an intermediate state 
m. This process produces characteristic logarithmic features in the tunneling 
spectrum. 


spin system has to be accounted for by the energy of the tunneling electron. The 
current flowing between tip and sample is then given by: 

-|-oo 

-ho 

-|-oo 

= /(e + ev)[l-/(e-e,j)]de, 

_ho 

(9) 

where we have dropped the additional summation that would account for current 
through the different sites in coupled spin systems, i. e. we restrict our calculation to 
single spin systems. 

The unitless tunnel barrier transmission coefficient Tq = \TtaJps\^ contains all 
experimental parameters that will determine the strength of the tunneling current. In 
particular, Tq strongly depends on the distance between tip and adsorbate, and on 
the spin independent local densities of states in tip and sample. For the moment we 
want to assume that Tq <C 1 so that the influence of the tip on the spin system is 
negligible and the time between consecutive tunneling events is much longer than the 
relaxation time of the spin system, e. g. via processes such as described by equation[8l 
Such conditions are usually easy to provide in STM experiments where the tip-sample 
separation can be adjusted to give tunneling currents in the pA range or below so that 
the tip can be seen as a local probe of the spin system that leaves the spins always in 
thermal equilibrium with the substrate (zero-current approximation). 

Under this assumption the average state occupation Pi{T) of the spin system is 















10 


only governed by the effective temperature T and follows the Boltzmann distribution, 


MT) 


exp[-ej {ksT)] 
Eiexp[-ei/(fcBT)]’ 


( 10 ) 


with ks as the Boltzmann constant. Furthermore, we assume a flat density of states in 
tip and sample so that the electron occupation is given by the Fermi-Dirac distribution 
/(e) = [1 + exp(e/fcBr)]-b 

While most STM experiments do not record the tunneling current directly, 
we are interested in its derivative with respect to the bias voltage. With the 
energy independent matrix elements the derivative of the current is easy to evaluate. 
Integrating the Fermi-Dirac distributions of equation |9] and calculating the derivative 
results in a temperature broadened step function [98]: 


0(a:) 


1 -I- (x — 1) exp(a;) 
exp(a:)^ 


( 11 ) 


with X = €/{kBT), so that equation [HI becomes: 


w 




27re^ 


rp2 

-^0 


'i'J 




0(eF - Cif). 


( 12 ) 


The main goal of this paper lies in deriving the transition matrix elements Mif to 
be able to solve the equationfT^ First, we start with second order processes, neglecting 
scattering that involves electrons which originate and end in the substrate bath, i. e. 
we concentrate on tunneling processes which consist of only one scattering event as 
depicted in figure ISb and described by the transport Hamiltonian of equations |S| and 

^ 1 . . 

= (</'/>V'/l 2^-S + = Mif + USif. (13) 

This matrix is energy independent and connects the initial states in the electron 
baths \Lpi) and spin system |'0i) with their final states \'ff) and \ipf)- It contains the 
spin-exchange scattering term Mif and a potential scattering term that has non-zero 
matrix-elements only when initial and final spin state are the same is the delta 
distribution). Computing the absolute square of results in three terms that 

contribute to the tunneling conductance: 


M 


( 1 ) 

if 


2 


+ \U\^d,f + 25?([/ X M,/)%, 


(14) 


where 5ft(H) denotes the real part of the matrix A. The first term consist of the 
operators 

(T • S — *^-t- ~b i^zSz^ (1^) 

which connect the initial and final state and accounts for spin-exchange processes in 
which angular momentum between the tunneling electron and the localized spin system 
can be exchanged. The second term accounts for potential scattering between the 
tunneling electron and the localized spin system and does not change the spin state. 
The third term results from the interference between potential and spin-exchange 
scattering and depends, as we will see, strongly on the angular momentum of the 
localized spin and is the origin of magneto-resistive tunneling (TS] |99] . 
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Due to the product-state of the total quantum-mechanical system, it is worth 
mentioning that the matrix elements have to be independently evaluated for the 
localized spin and for the tunneling electron. For an electron tunneling between a 
spin-average tip and sample the non-zero matrix elements are: 

(a) (t |s+| 4-) = -1-1, (c) (t Is^l t) = +1/2, (e) (t |/| t) = +1, 

(b) ais_it) = +i, (d) (iis,u) = -i/2, (f)u|ju) = +i, (16) 

where, for completeness, we have additionally included the probability amplitudes for 
interacting with the unity operator I. Since a = 2s, these matrix elements enable us 
to rewrite the spin-exchange scattering term in equation [TJ] for a spin-average tip and 
sample to [8l[99]: 


^ {ipf\S- \^i) 


1 

■^2 


{ipflS+l'ipi) + 


(17) 


Note, that mt and [TBh do not cancel out because the two different spin directions 
in the initial (tip or sample) and final (sample or tip) bath are incoherent ensemble 
states and therefore cannot interfere with each other. 


For a tip with spin polarization t] = 


pt—pi 

Pt+Pi 


along the z-axis, the transition matrix 


elements are tunnel-direction dependent because either the initial or the final electron 
bath has now a spin imbalance in the density of states [33] : 

1 + ry 




|M| 


s— 1 2 


2 

1 — ?7 


{ipflS- It/,) 

2 1 — 77 

+ 2 

(t/zI 5+ |7/i) 

2 

+ 

(V'/l IV'i) 


2 1 -hT? 

2 

(t/zI 5+ |7/i) 

2 

(t/zI \'ipi) 


(18) 


As an example figure 0] shows the spectrum for a single spin with S = 1/2. When 
applying the magnetic field, the conductance increases step-like at bias energies above 
the Zeeman-energy |eF| > g^B\B\- These steps are symmetrical around zero voltage 
and thermally smeared out by about 5 /cbT [^. When the electrons in the tip are 
spin-polarized the step heights at positive and negative bias are different. For a 
paramagnetic tip with a spin polarization in direction of the applied field (ry > 0 ) 
the conductance step at positive bias decreases, while it increases by the same amount 
at negative bias (Figure |4 )d). 

With a spin-polarized tip the interference between the potential and Kondo- 
scattering shows a very particular effect [H 1^ . Depending on the sign of U the 
zero-bias conductance is either increased or reduced and can become even zero at 
77 = 1 and U = —(Figure |31:-d). This magneto-resistive elastic tunneling [51] 
arises from the third term in equation 1141 and scales with the expectation value {S^) 
of the impurity. At non-equilibrium the strength and sign of this term can change; 
this allows to read-out of the z-projection of the spin without exciting it, as has been 
shown in spin-pumping (see section El) [13199] and pump-probe experiments |14j . 


2.2. Expansion to third order 

We now turn our attention to the next higher order of interaction. We want to consider 
processes in which, during the transport of an electron from tip to sample (or vice 
versa), the spin system additionally interacts with an electron from the sample bath 
(Figure Et)- These processes involve an intermediate state (m) and are expressed in 
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(a) (b) 




Figure 4. Tunneling spectra for a 5 = 1/2 system at a temperature of T = 1 K 
and g = 2 evaluated in second order perturbation theory. (a) At applied 
magnetic field temperature broadened symmetric steps are visible. The different 
contributions to the conductance are labeled. Gray: conductance due to potential 
scattering {\U\ = 0.25), green: due to spin-spin scattering without changing the 
localized spin, yellow: spin-spin scattering with changing the spin, (b)—(d) A 
spin polarized tip produces asymmetric steps. Here, the zero voltage conductance 
depends strongly on the potential scattering U and the tip polarization. 


second order Born approximation as: 


M 


( 2 ) 

-i-f 


J Ps 





Cm — e/ 


(19) 


Here, the tilde on the matrix M tags scattering processes between the localized spin 
and sample electrons only. The order of the two different electron-spin interactions can 
be exchanged. Thus, we have also to account for processes in which first an electron 
that originates and ends in the sample scatters the system into the intermediate state 
and then the tunneling electron interacts with the system scattering it into its final 
state (right fraction in equation [19]). Here, it is of fundamental importance to have in 
mind that, contrary to the initial and final state of the total system, the intermediate 
state I'lpm) can be virtual, i. e., must not necessarily obey angular momentum and 
energy conservation. The integro-summation symbol in equation [JHI indicates that 
we perform both, a summation over all possible discrete intermediate states in the 
local spin system \ tjj) and an integration over the continuous states of the intermediate 
electron states \ip) in the substrate. For the integration, we consider states in an energy 
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range of zLloq around Ep as possible scatterer leading to the following characteristic 
function for electron-like processes [HZl 1100 ) : 

+ OD +W0 

Feie, T) = - I de" J de' ^ /'{e” - e, T), (20) 

— OO —tOQ 

and for hole-like processes: 

+ 00 

FH{e,T) = - J de" J de'{e" - e,T), (21) 

— OO —loq 

respectively. Here, the Fermi-Dirac distributions in the numerator ensure an 
unoccupied (occupied) final state and the integration over the derivative f accounts 
for the temperature broadening during the total process (Figure [SK). The switching 
from virtual to real processes at e = 0 leads to a logarithmic singularity that is 
broadened by the temperature. Thus, as long as the energy of the intermediate 
state and the temperature kpT are small compared to the integration bandwidth wq, 
the equations [20] and [21] (with a change of sign) can be rewritten as: 

+ 00 

F{eV -e^)^- [ de' In (^)- (^2) 

J \eF-em + Ho/ 

— OO 

Here, 0'(e,T) = dQ{e,T)/de is the derivative of the temperature broadened step 
function ('equation 111)1 and a small Fq accounts for additional non-thermal (lifetime) 
broadening. Figure [5] shows the energy and temperature dependence of F which 
has first been experimentally observed on a radical molecule adsorbed on a Au(lll) 
surface (see figure [T^ and f) [Hj. Compared to the Lorentzian or Frota functions 
|101) normally used for describing Kondo resonances in STM measurements [12], the 
behavior of this function is quite different, even though - if analyzed only in a small 
energy range of a few ksT around zero -, it has a similar shape [22] . The function 
shows a relatively flat top whose width is determined only by the temperature. For 
energies kpT < |e| < wq the amplitude decays logarithmically as — ln(|e/wo|) while 
the maximum peak height is proportional to —ln{kBT/uJo). Due to the restriction 
to states in an energy interval ±wo the function remains analytical even at e —>■ ±oo 
where F approaches zero. Note, that the precise value of the cut-off energy loq is not 
critical and mainly changes only the background offset. If not otherwise noted we use 
Wo = 20 meV throughout the paper. 

When calculating the conductance we now have to consider both processes 
depicted in Figure [3 ]d and c and thus have to replace A4i with Mi + M 2 Fv equation 
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Figure 5. (a) Scheme of the third order processes. Full (dashed) vertical arrows 
illustrate real (virtual) transitions between intermediate and final state (m —)• /) 
of electron-like (full circles) and hole-like (empty circles) carriers, (b) and (c) 
The temperature broadened logarithmic function F(eV — €jn-,T) at different 
temperatures T with cjq = 200 meV and Fq = 5 p.eV. (b) linear, (c) linear-log 
plot. 


[12] leading to: 




M 


( 1 ) 

if 


Jps F{erm) + 

m 


JpsU ^ F{emi) + 

m 


O {Jpsf 


(23a) 


(23b) 


The evaluation of the matrix elements up to third order yields two new terms due to the 
interference between the processes described by and A42 which are absent in the 
second order perturbation calculation. The term l23^ was first identified by Jun Kondo 
[34] and can lead to temperature broadened logarithmic features in the conductance at 
the energy of intermediate states. For a non-vanishing potential scattering amplitude 
(U ^ 0) the term l23bl will. in addition, produce a bias-asymmetry in the differential 
conductance even without spin-polarized electrodes. 

We start with the evaluation of the Kondo-like processes, that are described by the 
term l23^ using a spin S' = 1/2 system with only two states as an example. We assume 
that in thermal equilibrium only the ground state |t) = IV'i) is occupied, i. e. that 
the Zeeman-splitting induced by the external magnetic field is large compared to the 
thermal energy: gpLB\B\ ^ ksT. The Feynman diagrams of figure ]6] depict all possible 
interaction processes. First, we review the second order processes (Figure]®^): A spin- 
up electron that tunnels from tip to the sample cannot flip the spin (|'0i) —>■ IV'i))) 
while a spin-down electron can either scatter with exchange of angular momentum 
(I'i/’i) —>■ |V^ 2 )) leaving the system in the state H) = \ip 2 ) or without exchange 
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(a) 






Q{Ei-E^+eV)XF(Ei-E^+eV) Q{E-Ef+eV)xF{E^-E.+eV) 


Figure 6. Interaction diagrams of order two (a) and three (b) for an electron 
tunneling from tip to sample into a two level 5 = 1/2 spin system. The large 
(small) spheres depict the state of localized spin (interaction electron) and the 
color their spin directions. Schematic spectra show their contributions to the 
conductance at positive bias. The numbers label the processes with the state 
order of the localized spin-system. An appended ’R’ label processes in which the 
scattering into the intermediate state is performed before the tunneling electron 
interacts (exchange diagrams). Note that the time order of the processes influences 
crucially the conductance spectra as schematically displayed in the small graphs 
(vertical line is Ep; the * means multiplication). 


(I^i) —^ IV'i))- To third order, there are a total of six diagrams to be accounted 
for, which we can label by the states occupied in the initial (i), intermediate (m), and 
final (/) state. Here, we mark the exchange diagrams by appending a ’R’. This results 
in the processes (112), (121), (122), and the reversed order scattering events (112R), 
(121R), (122R). 

To evaluate the conductance due to these processes we calculate the spin-flip 
matrix elements of equation [15] for the electrons and the localized spin-system. For 
electrodes without spin-polarization this results in conductances of: 




47re^ 

h 


'^0^ -5; IV'/) (V'/l IV'm) (V'n 


j,k,l = 


EL 

dV 




xF(eV- 

, JPs 


47re^ 


-T, 


4z 


,T)x0(eR-e,/,T), 

£jkl^ (V’jI IV'/) (V'/l IV'm) (V'n 


E 

j,k,l = 

{x,y,z} 


Sj \tpi) 


Sj \tpi) 


(24) 


xF(eV-e™„T) x 0(eR-e,/,r). 


( 25 ) 
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(a) (b) 



Figure 7. Tunneling spectra for a 5 = 1/2 system (a) The third order 
contributions for non spin-polarized tip and sample with ksT = 0.1 x and 

cjQ = 20 X gfisB. (b) Total conductance for Jps = —0.05, U = 0, T = 1 K, and 
g = 2 at different magnetic fields, (c) Additional contributions to the conductance 
in third order due to a non-zero potential scattering term, (d) Total conductance 
as in (b) but with U = 0.25. Dashed lines in (b) and (d) show the second order 
conductance only at B = 10 T. 


Here, equation accounts for the direct diagrams (normal order) and [551 for the 
exchange diagrams (reversed order), respectively. In these equations ejki is the usual 
Levi-Civita tensor of rank three, which is 1 (-1) if {jkl} is an even (odd) permutation of 
{xyz}, and zero otherwise [77]. The step-function 0 ensures that energy conservation 
at the final state of the scattering process is obeyed, while the function F has its peak 
at the intermediate state energy and is mirrored at Ep for the reversed tunneling 
process. Any spin-polarization in tip or sample changes the transition matrix elements 
for the interacting electrons and equations [53] and [55] would become more complex (see 
[Appendix Aj ). 

Comparing the different contributions to the conductance (figure [Hd and [7^) it is 
remarkable that all third order scattering events start with a spin-down electron except 
the process 121i?. Here, a spin-up electron tunneling from tip to sample non-trivially 
interacts with the system because a substrate electron has flipped the localized spin 
into the intermediate spin-down state before the interaction takes place. 

Summing up the contributions to the conductance due to all second and third 
order scattering events results in spectra similar to those exemplary given in figure 
[7b (c. f. figure [lb and f [22]). At zero field the differential conductance shows a peak 
that splits with increasing magnetic field. While in this calculation we assume to be 
in the weak-coupling Kondo limit and thus neglect any correlation energy due to the 
formation of a Kondo singlet, the peak splits as soon as the Zeeman energy overcomes 
the thermal energy. Note, that the resulting split-peak at small fields can lead to the 
deduction of an erroneously high g-factor when just evaluating the peak positions due 
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to the superposition of peak and step-structures. 

If we consider now in addition that the scattering process between tip and sample 
or vice versa can also occur via the potential interaction Iterm [23bl) . we observe an 
asymmetric line-shape as soon as the two eigenstates are no longer degenerate and 
(Sz) ^ 0 (figure [3: and d). This direction-dependent asymmetry cannot originate 
from the scattering matrix elements that involve only the localized spin site (the 
order of excitations shall be the same) but derives from the matrix elements involving 
the interacting electrons. Physically, the reason lies in the asymmetry of the tunnel- 
junction for which we assume that only the sample is coupled to the spin-system and 
thus neglect any intermediate scattering process that originate and end in the tip. 

As an example, we examine in detail the process (121) (see figure [6 )d) in both 
tunneling directions. For a current to flow from tip to sample via this process, first a 
14,) electron originating from the tip is scattered into a |t) state in the sample exciting 
the spin system from {ipi) —>■ \ijj 2 )- Second, a |t) electron in the sample is scattered 
into 1),) bringing the spin system back to \ipi)- Third, the sample jj,) is scattered back 
into the tip as jj,) electron. The last step of this process can either take place via 
the (7z or, in the case of potential scattering, via the d/ term, respectively. For an 
electron tunneling in the reverse direction we first have a jj,) from the sample being 
scattered into an electron |t) state of the tip simultaneously exciting the spin system 
from l^i) —>■ \ip 2 )- Then, the second process flips a |t) sample electron into the ||) 
hole from the first process. Finally, a |t) tip electron traverses the junction and fills 
the It) hole in the sample. Comparing both tunneling directions we see that it is the 
third process which differs in the initial and final state, i. e. the matrix elements are 
either (t \az\ i) and (t |d/| t) or (f \az\ t) and (f |(T/| f)- Rearranging the matrix 
elements so that all processes become electron-like, the prefactors for calculating the 
tunneling conductance due to these four processes are for the two without potential 
scattering: 


(121)*^^ : + Jp, (i‘ |dz| i^) (i^ |d_| f) (f Id+I i*) > 0 

- + + 

(121)^^* : - Jps (f |d+| i^) (r |<T_| r) (r |dz| f) > 0, 

+ + + 


and with potential scattering: 

( 12114 )*^" : 

(121u)"^* : 


+ Jp.t/U* \di\ i^) (r |d_| r) (f |d+| i‘) < 0 

+ + + 

- jpsUit |d+| r) (r |d_| f) (f |dz| f) > 0, 

-I- -I- -I- 


(26) 


(27) 


where we assumed Jps < 0 and U > 0. Note, that the preceding sign change at 
the tunneling direction from sample to tip (s —> t) is due to the rearrangement 
of the interaction order together with the switching from hole-like to electron-like 
scattering. The contribution to the conductance from the processes in which only 
Kondo-like spin-spin interactions take place is positive for both tunneling directions, 
while the conductance for processes that include potential scattering changes its sign 
when inverting the tunneling direction. 
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3. Some single spin examples 

Up to here we have revised the relevant formalism to calculate the conductance in 
the third order of the matrix elements quite closely following the work established 
by Appelbaum, Anderson, and Kondo [MJ [351 [Ml 122] ■ The S = 1/2 system we used 
for illustration contained only two states and was not influenced by any magnetic 
anisotropy or near-by spin systems. Recently, efforts have been made to expand this 
perturbative model to higher spin systems, which also include magnetic anisotropy and 
couplings to neighboring spins IHIHSl [HI [SH, but the importance of the potential 
scattering was not taken into account up to now. In the following, the power of this 
easily accessible model will be used to evaluate and describe the spectral features on 
more complex systems. 

3.1. The smallest (S =1) high-spin system 

Regarding the spectra of the S' = 1/2 system (Figure [7|), one could get the impression 
that the coupling to the sample always results in some superimposed peak-like 
structures as soon as a spin flip transition is possible, and which scales with the 
substrate coupling strength —Jps. In this section we will show that the situation is 
more complex and depend not only on the transition matrix elements but also on the 
state energies. 

To expand the complexity, we turn now to a magnetic system with S' = 1 in which 
axial and transverse magnetic anisotropy have broken the zero-field degeneracy of the 
three eigenstates. Assuming easy-axis anisotropy {D < 0) the energetically lowest 
eigenstates have weights only for = ±1. The additional transverse anisotropy splits 
the remaining two-fold degeneracy forming an antisymmetric ground and a symmetric 
first excited state (Figure [HK). Thus, in the basis of the magnetic easy axis the three 
eigenstates can be written as: = -^(1 + 1) “ l~l)); IV^2) = "^(l + l) + 

IV's) = |0)- Such situations have been found for example for Fe(II)-phthalocyanine 
molecules adsorbed on the (2 x 1) oxygen terminated Cu(llO) surface [65]! for Fe atoms 
adsorbed on InSb(llO) [53], or for CoH complexes adsorbed on the /i-BN/Rh(lll) 
surface m- 

By taking only second order spin-flip processes into account it is easy to calculate 
(using equation [T7|) that the transition probabilities from the groundstate to the two 
excited states at R = 0 are equal, leading to a symmetric double step structure in 
the spectrum (Figure jH]). Expanding the calculation to third order, only the processes 
which involve all states, i. e. (123) and (132), have non-zero matrix elements when 
the potential scattering [7 = 0. Processes like (112), (121), or /131) cannot be 
interlinked using any combination of the operators 5”+, S-, and Sz and are thus 
forbidden. Nevertheless, additional potential scattering enables transitions involving 
the operators S+, S-, and I or Sz, Sz, and I which makes the aforementioned processes 
possible resulting in an asymmetry of the spectrum with respect to tunneling direction. 

It is remarkable that the remaining third order processes at U = 0 change the 
spectrum in a quite different fashion, even though they have the same strength. 
Process (123) produces its peaks at eV = ±e 2 , but due to energy conservation it 
contributes to the spectrum only at |eU| > ea, which efficiently cuts off the peak. In 
contrast, the peak at higher energy due to process (132) is not as strongly cut off 
(Figure [Hb). Thus, the full spectrum shows a peak-like increase of the conductance at 
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Figure 8. Tunneling spectra of a spin S = 1 system, (a) The anisotropy lifts 
all degeneracies [D = —5 meV, E = 1 meV) at zero field. The plot shows the 
state energy over the expectation value (m^) = tr (^Sz {ipi) (V'il) (blue: E = 0, 
black: E ^ 0) (b) The significant third order processes (T = 1 K, iJ = 0 T). 
The dot dashed line show the additional contribution for a U = 0.125. (c) The 
total spectrum with Jps = 0 (dashed red line), —0.1 (full black line), and with 
additional U = 0.125 (dotted dashed blue line). 


the energy of the second step but not at the energetically lower first transition step. 
Furthermore, the conductance has a curved form for tunneling voltages between e2 
and ea (Figure (5}:). Note, the overall general behavior would not alter when changing 
from easy-axis to easy-plane anisotropy {D > 0) as long as all degeneracies are lifted. 

3.2. Single Mn and Fe atoms on CU 2 N 

After the, with only three eigenstates, rather easy S' = 1 example, we now apply the 
model to the experimentally and theoretically intensively studied single 3d transition 
metal atoms adsorbed on a monolayer of CU2N on Cu(lOO). When these atoms are 
placed on top of a Cu site they form strong covalent bonds with the neighboring 
N atoms [H]. This highly anisotropic adsorption geometry leads to three distinct 
symmetry axes that are perpendicular to each other: The direction out-of-surface and 
two in-surface directions along the Cu-N bonds and perpendicular to it, along the so 
called vacancy rows (Figure IHt- inset). 

Single Fe atoms adsorbed on this surface have been found to be in the S = 2 
state with a magnetic easy-axis along the N rows (z-direction) and a magnetically 
hard-axis along the vacancy row (a;-direction). In this coordinate system anisotropy 
values oi D = —1.55 mV and E = 0.31 mV, and a gyromagnetic factor of 5 = 2.11 
described the experimental data well using a spin Hamiltonian like equation [S] and 
a second order tunneling model [H [75l 1102] . The Hamiltonian has as solution five 
non-degenerate eigenstates and, due to Z? < 0, favors, at zero field, ground states with 
weights at high values. Similar to the 5 = 1 system the transverse anisotropy 
breaks the degeneracies leading to a symmetric and antisymmetric solution with the 
main weights at |±2) as ground and first excited state and weights in |±1) for the 
second and third excited state (Figure IHlo). To visualize the anisotropy we plot the 
total energy necessary to rotate the ground state into arbitrary directions showing the 
favored easy axis (z) and the unfavored hard axis (a;) (FigureEb) [89]. In second order, 
spin-flip scattering is allowed between the groundstate and the three lowest excited 



















dl/dV (arb. units) dl/dV (arb. units) 


20 



bias V (mV) 


bias V (mV) 


Figure 9. Comparison of experimental and calculated spectra on Fe and Mn 
atoms adsorbed on Cu 2 N. (a) Experimental data from ref. [8] measured on a 
single Fe atom at increasing field Bz along the easy axis and at a temperature 
T = 550 mK (red circles). The simulations (black lines) for all plots are obtained 
with one set of parameters: g = 2.11, D = —1.57 meV, E = 0.31 meV, 
Jps = —0.087, U = 0.35, and Tgff = 740 mK. Additionally, a constant offset 
of 20% of the total conductance had to been added (black shaded area). For 
Bz = 0 the second (green shaded area) and third order (orange hatched area) 
contributions to the conductance are indicated. The spectra at field are vertically 
shifted for better visibility. The inset shows the adsorption site of the 3d atoms 
(black circle) on the Cu 2 N (Cu yellow, N blue circle) (b) and (c) Schematic state 
diagrams and visualizations of the magnetic anisotropy (in meV) for Fe (b) and 
Mn (c) |89| . (d) Experimental data of two different Mn atoms at Bz = 0 and 
Bz = 7 T (colored circles). The fits (full lines) for the Bz = 0(7)T data results 
in Jps = -0.029(-0.0091), U = 1.35(1.28), D = -51(-39) //eV, g = 1.9, and 
Teff = 790(930) mK. The dashed line simulates for the zero field data the absence 
of any anisotropy. The dotted line simulates with the 7 T parameters the absence 
of a magnetic field, (e) The 7 T atom as in (d) probed with tips of different spin 
polarizations 77 * (Data from ref. [ 99 ]). 
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states but a transition to the highest state is forbidden because this would require an 
exchange of Am = ±2. 

Experimental dl/dV measurements on this system show, in addition to the 
conductance steps, peak-like structures at the second and third step but not at the 
lowest one (Figure |9^). Additionally, they show an asymmetry between positive and 
negative bias. To rationalize these observations we can follow a similar argumentation 
as in the S' = 1 case: In third order, transitions like (121) are not possible without 
additional potential interaction and processes like (123) or (124) are strongly cut 
off due to the high energy difference between €2 and £3 or £4. In contrast, the 
processes (132) and (142) scale with Jps leading to the peak features in the differential 
conductance. The additional asymmetry hints at a non-negligible potential scattering. 
As the computed curves in figure [5^ reveal, our model almost perfectly fits the 
magnetic field data without any adaption of the parameters^ We mention that 
similar good agreement between experimental data and computation can be reached 
in the other two magnetic field directions (not shown). The coupling strength in these 
simulations is Jps = —0.087, close to the —0.1 found in a similar perturbative approach 
|84] . A potential scattering term oiU = 0.35 is necessary to reproduce the asymmetry. 
This value is significantly smaller than the U ~ 0.75 found in experiments where the 
magneto-resistive elastic tunneling was probed [99) . Part of this discrepancy can be 
understood by an additional conductance term that does not coherently interact with 
the spin-system and which would lead to an overestimation of U in magneto-resistive 
measurements. Indeed we need a constant conductance offset of about 20%, which is 
added to the calculated conductance to reproduce the spectra. 

Switching from an integer to a half-integer spin system we now discuss individual 
Mn atoms on CU2N, which have a spin of S' = 5/2 and only a small easy-axis anisotropy 
of H ~ —40 peV along the out-of-plane direction and a negligible transverse anisotropy 
[ail]. The easy-axis anisotropy prohibits the immediate formation of a Kondo state 
due to a Kramer’s degenerate ground state doublet with ruz = ±5/2 (Figure |9t). At 
zero held a typical spectrum shows only one step, which belongs to the transition 
between the ±5/2 and the ±3/2 states that have superimposed asymmetric peak 
structures (Figure [5)1). The ht to the model yields Jps = —0.029 and resembles a 
5=1/2 split-Kondo peak at small magnetic Helds (see HgurejTlD). A different Mn atom 
investigated at i?2 = 7 T shows a signihcantly reduced Jps = —0.0091. Interestingly, 
we Hnd for both atoms a potential scattering value of [/ ~ ^5, which allows one to 
describe the spectra without the need of any additional conductance offset. This high 
U value that is the origin of the bias asymmetry has been independently found in 
spin-pumping experiments m ( see section [5]) and by measuring the magneto-resistive 
elastic tunneling contribution [99) . The extraordinary agreement between model and 
experiment can be seen in measurements using different spin-polarized tips on the same 
atom (Figure [5^). Here, the strong inhuence of the tip-polarization on the inelastic 
conductance at bias voltages |F| < 1 mV is evident while the differential conductance 
at V >1.5 mV stays, for all tips, constant. 

§ We note that all experimental data in this manuscript are unprocessed except for a slight slope 
adjustment of < 2% which corresponds to the compensation of an unavoidable small drift in the tip 
sample separation of less than 0.5 pm during the measurement time. The effective temperature is 
always larger than the base temperature of the experiments due to additional broadening effects as 
for example modulation voltages, electromagnetic noise, or additional lifetime broadening. 
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Figure 10. The tunneling spectra of Co atoms {S = 3/2) on CU 2 N. (a) Schematic 
state diagram and visualizations of the magnetic anisotropy (in meV). At B = 0 
the states \ip)i and 1-0)2 are degenerated and differ by Arriz = ibl. (b) and (c) 
Experimental data from reference m of two different Co atoms at = 0, 4,6 T 
(colored circles, top to down, shifted for clarity). The best fits (full lines) results 
in Tefj = 4 K, D = 3.3 meV, E = 0.7 meV, g = 2, Jpg = —0.11 for (b) and 
Teff = 3.8 K, D = 2.7 meV, E = 0.5 meV, g = 2, Jpg = —0.25 for (c), respectively. 


4. The limit of the perturbative approach: The Kondo system Co on 
CusN 

When a half-integer spin with 5 > 1/2 has easy-plane anisotropy D > 0, its ground 
state at zero field is a doublet with its main weights in = ±1/2. This enables an 
effective scattering with the substrate electrons and leads at low enough temperatures 
to the formation of a Kondo state. Experimentally this has been observed for Co atoms 
on CU 2 N 0 nn m nil [m which have been found to have the total spin S' = 3/2 
and which enter the correlated Kondo state with a characteristic Kondo temperature 
of Tx = 2.6 K in experiments performed on small patches of CU 2 N at temperatures 
down to T = 550 mK 0. Apart from D > 0 the system also has a small in-plane 
anisotropy {E ^ 0) which creates an easy axis (x) along the nitrogen row and a hard 
axis (z) along the vacancy rows (Figure fTUb). 

Interestingly, in this system the coupling to the substrate Jps changes with the 
position of the Co atom on larger CU 2 N patches, concomitant with a change in the 
anisotropy energies which separates the |±l/2) states from the |±3/2) ones [T6l|94j. 
For us this allows the study of the transition from the weak coupling to the strong 
coupling regime. In the case where the Co atom is relatively weakly coupled to the 
substrate {Jps ~ —0.1), the model can be consistently fitted to the experimental data 
even at different field strengths along Bx (Figure [TObl. We observe a zero-energy 
peak that splits at applied fields similar to the S' = 1/2 system of figure [7li. But while 
for S = 1/2 the field direction does not play a role, here the peak splitting depends 
strongly on the direction due to the magnetic anisotropy 0 . For this high-spin system 
we furthermore observe inelastic steps due to the transition to energetically higher 
excited states which are located at \eV\ = 2|Z1| for B — 0 and whose additional peak 
structure is well described within the model. 

For Co atoms adsorbed closer to the edges of the CU 2 N patches, the coupling 
to the substrate increases and the fit to the model worsens significantly (Figure [TOt). 
While the experimental data measured for non-zero fields are reasonable well described 
with Jps ~ —0.25, at B = 0 the experimentally detected peak at Ep is stronger than 
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the peak created by the model. Additionally, the experimental peak-width appears to 
be smaller than the temperature broadened logarithmic function, which was similarly 
observed for a radical molecule with S' = 1/2 for temperatures presumably close to 

Furthermore, we observe that the calculated spectrum no longer well describes 
the steps at 2|Z1|. At this energy the third order contributions are less pronounced in 
the experimental data, indicating that we reach the limit of the perturbative approach. 

The full description of a spin system in the strong coupling regime requires 
complex theoretical methods like numerical renormalization group theory |1031l36lll04j 
which are beyond the scope of this paper. Nevertheless, we can discuss some of 
the physical consequences within the outlined perturbative model. In contrast to 
the examples discussed in section [3l the two lowest degenerate ground states of the 
C 0 /CU 2 N system have weights in states that are separated by Am = ±1. Thus, at 
zero field, electrons from the substrate can efficiently flip between these two states. 
The computation of the transition rate between the two states \’4>i) and |'0/) of the spin 
system that have the energies and e/ is := hf/e, i. e. the transition probability 
per time, is analogous to the calculation of the current (equation given by: 

cy nOO 

= ^ = J de\Ml/^\^f{e,T)[l-f{e-e,f,T)]. (28) 

Here we have set the bias to eF = 0 and replaced the coupling constant Tq between 
adsorbate and tip with Jps, the coupling between adsorbate and substrate. When we 
now consider only processes up to second order, the matrix Mif is independent of the 
energy and the integral can be evaluated to 

de/(e,T)[l-/(e-e,/,T)] =--. (29) 

Thus, the scattering rate between the two degenerate states (i. e. ei 2 = 0) is directly 
proportional to the temperature (Figure [TT^): 

^(2)s^s ^ (30) 

These scattering processes, which we have discussed in section [5] and displayed 
in figure |6^, will tend to change the spin polarization of the electronic states in the 
sample near the adsorbate to be correlated with the localized spin. Nevertheless, this 
local correlation will be quickly destroyed by decoherent scattering processes with 
the remaining electron bath, which we can assume to be large and dissipative. This 
decoherence rate is also proportional to the temperature, Tdecoh oc ksT |105| . but 
usually stronger so that no highly correlated state can form. This picture changes when 
we additionally consider third order scattering processes which yield, using equation 
Uni the probability: 

A poo 

= -riJp.f / de 

X [F{erm - e, T) + F(e,™ - e, T)] /(e, r)[l - /(e - e./, T)]. (31) 

Due to the growing intensity of F(e = 0) at reduced temperatures (Figure O, for 
temperatures T —0 the scattering F^^^ decreases significantly more slowly than F^^^ 
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(a) 


(c) 




(b) 


E (meV) 


(d) 




e (meV) eV (meV) 


Figure 11. Correlations induced by substrate electrons, (a) In second order the 
scattering probability of a substrate electron with energy e is given by the overlap 
of the electron and hole-like Fermi-Dirac distributions (area underneath the curve) 
and is for degenerated ground states directly proportional to the temperature. 
Third order scattering (b) gets logarithmically stronger than the second order 
processes with decreasing temperature, (c) Schematics of the entangled state at 
temperatures below the characteristic Kondo temperature. The scattering leads 
to a correlation between the spin state of the localized system and the sample 
electrons in the vicinity, (d) Spectrum of a Co adatom on CU 2 N in the Kondo 
regime {B = 0 T, T = 0.55 K, data from ref. [S])- The experimental dlldV curve 
(red dots) is well described by the sum (full line) of a Frota peak centered at 
Ep (dotted dashed line) and a broadened step-function at higher energy (dashed 
line). 


(see figure [Hb) so that their ratio steadily increases: 


r(3) 

r(2) 



(32) 


In contrast to that, the decoherent scattering rate with the bath lacks localized 

scattering centers and therefore has no significant third order contributions. Equation 
[5^ leads to a characteristic temperature, the Kondo temperature Tr-, where and 
higher order scattering terms become the dominant processes and perturbation theory 
breaks down [34l 1106111071 136] : 


Tk 


Wo f I \ 


(33) 


Below this temperature the assumption used up to now, i. e. that the electronic 
states in the sample are not influenced by the presence of the spin system, no longer 
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applies. Using exact methods like the modified Bethe ansatz |108l I109j or numerical 
renormalization group theory [104] reveals that the sample electrons in the vicinity 
rather form an entangled state with the impurity, i. e. 




total 




(34) 


as illustrated in figure ITTb . This combined state is quite complicated because the 
electronic states in the sample are continuous in energy and extend spatially, and 
therefore strongly alter the excitation spectrum of the adsorbate [11011111] . 

Figure ITTH shows the impact of the formation of the correlated state on the 
experimentally detected spectrum of the C 0 /CU 2 N system. At temperatures T Tk, 
the peak at the Fermi energy can no longer be well reproduced by a temperature 
broadened logarithmic function (which in any case must diverge for T —>■ 0). The 
peak is much better described by an asymmetric Lorentzian, i. e. a Fano line-shape 
[47] , or the Frota function [1011111211113] which has a finite amplitude and a half¬ 
width of A = ksTK, corresponding to the correlation energy of the Kondo state. 
Superimposed on this zero-energy peak we detect conductance steps at voltages that 
enable scattering of the spin system into the = ±3/2 state. These steps are well 
described using only a second order perturbation spin-flip model |9] and do not show 
any third order logarithmic contributions. This means that the probability of the 
third order scattering channels must be closed due to the ground-state correlation 
between the localized spin and the substrate electrons. Here we see that the simple 
perturbative approach of the model is no longer valid and fails to capture the full 
physics in this strongly coupled Kondo regime. Nevertheless, the model still gives 
valuable information because it enables us to identify the conditions under which 
correlations can evolve and is due to its computational simplicity also applicable for 
larger and more complex spin systems where the calculation of the exact solution 
might not be feasible or very time consuming. 


5. Spin dynamics and rate equations 

In the preceding section we discussed the appearance of correlations due to higher 
order scattering between the substrate electrons and the localized spin system. Now 
we want to evaluate how the tunneling of electrons between the biased tip and the 
sample influences the state populations and the observable spectroscopic features in 
local conductance measurements. This is equivalent to dropping the zero-current 
approximation we have assumed until now. We will see that the change of time- 
averaged system properties results in characteristic fingerprints in the dl/dV signal 
which can become crucial for a deeper understanding of the system’s response and 
dynamics. 

To describe the transition rate between the localized spin states {ipi) and lijjf) 
due to the tunneling current flowing between tip and sample we can again use 
equation El and the relation r)j^®(eU) := /)/^®(eU)/e, which leads, in first order Born 
approximation, to: 

O-TT /*00 

= —T^ /(e- eU,r)[l - /(e- e,^,r)]. (35) 
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The integral can be solved using equation [55] resulting in: 


-pi (2)t— 


ieV) 



t—^s |2 

if I 


exp 



- 1 


(36) 


As long as eV — Cif ^ ksT, equation |36| can be further simplified to an equation that 
is linear in the energy difference: 


rif = {2^/h)Tl\M^^^\\eV - e./). (37) 


This linear dependence of the scattering rate on the energy difference is equivalent 
to the assumptions that in second order perturbation the matrix elements and the 
coupling constant Tq between tip and sample are energy independent. Under these 
assumptions the rate is given by the energy window between the available energetically 
hot electrons and the energy needed to change the localized state m- However, one 
should keep in mind that at large energy differences the intrinsic variation of the local 
density of electronic states in tip and sample p(e)*’® might alter the results. 

Furthermore, we additionally include third order contributions to the tunneling 
transition rates. This can be archived by integrating the interference terms between 
M.\ and Ad 2 of equation [23l and leads to an additional contribution to the scattering 
rate of: 

\ PCO 

r!f ^^(eU) = -rT^Jps / de ^ /(e - eU, T)[l - /(e - e.f,T)] x 

^ ^-00 ^ 


T) + 
T) + 
^mi T) + 
U^{If,M^fM,^)F{e,^ -e,T) . 


(38) 


For tunneling currents flowing in the opposite direction, i. e. for electrons which are 
scattered at the spin system when tunneling from sample to tip, equations 1361 and 1381 
can be adapted straightforwardly. 

The total transition rates F‘j^® = and F®”^* = 

have to be evaluated for all possible initial and final localized spin states and, together 
with the dissipative substrate to substrate scattering rates discussed in section [4] 
(equations [28] and |3T]) , can be brought together to the characteristic master equation 
for the state populations: 

= J2p,{t){T]r + rr‘ + + rr* + 

j¥^i j¥^i 

= Tij{t)pj{t). (39) 


This set of first order differential equations describes the time-development of an initial 
state of the spin system, which might be a superposition state |V'(0) = J2iPii^) IV'j)) 
at the time t under the influence of all possible, bias dependent, scattering events. The 
first summation in equation|39|accounts for all the probabilities to scatter into the f-th 
eigenstate from any other eigenstate. These probabilities have to be weighted with the 
time-dependent state probability pj of the originating state. The second summation 
accounts for all scattering events which reduces the f-th state probability by scattering 






27 


it into other eigenstates. It is important to note, that we use three crucial premisses 
in this approach: 


(i) We assume that the future of the localized quantum state depends only on 
the actual state and not on any interactions which happened in the past. 
This Markovian limit means that we assume that the timescales at which this 
approximation breaks down, in particular where higher than two-time correlations 
are dominant, can be neglected [114] . 

(ii) We do not account for any phase coherences between the eigenstates of the spin 
system. This means that the phase coherence time is short compared to the 
state’s lifetime r and the time between successive tunneling events. Under these 
assumptions all off-diagonal elements of the density matrix in the rotating frame 
with the eigenstates ip' are zero and thus the density matrix can be described as 
a state mixture X = J2iPi IV’i) (V’il- 

(iii) Because we treat the total system as the product state between the continuous 
electronic states in tip and sample and the discrete spin states, a correlated 
state between substrate electrons and the localized spin as discussed in section |4] 
(equation 1341) can not be developed directly within the limitations of our model. 


In the following, we further want to limit our evaluation to steady-state conditions, 
which means that the change of any external parameter, like the tunneling bias voltage, 
occurs adiabatically slowly, much more slowly than any relaxation times in the system. 
The steady-state condition is reached, when all pi are static, i. e. —)> oo) = 0. 

This is equivalent to finding of the algebraic kernel of the rate matrix T, which 
additionally has to be normalized to account for the conservation of probabilities 


Y.^Pi = 1 : 




kerT(eU) 

||kerT(eU)||- 


(40) 


From these p®*“*(eU) the current I{eV) and the differential conductance dl(eV) /d{eV) 
can be calculated: 


/(eU) = e^pf“‘(eU) (r‘;"^(eU) - q;"‘(eU)) . (41) 

ij 

To show the influence of a non-zero tip-sample coupling strength Tq on the 
tunneling spectra we are returning to the spin 5 = 1 example from section 13.11 
In figure IT^ and c we show the development of the dl/dV spectra a.t B = 0 and 
= 10 T for different coupling strengths between sample and tip, respectively. The 
coupling constants correspond to a tunneling resistance between Rt ~ 130 Mil and 
Ki 2.6 Mil, which is equivalent to a tunneling current of / ~ 150 pA - 8 nA at a bias 
of U = 20 mV. Note, that these are quite typical parameters for STM experiments. 

Without magnetic field and at small couplings we observe a spectrum which does 
not differ significantly from the one calculated in the zero-current approximation. The 
average occupation of the two excited states increases only slowly for voltages above 
the threshold for these transitions and reaches at U = ±10 mV only about 4% for p 2 
and less than 0.4% for p^ (dashed lines in figure IT^l . This situation changes when 
the tip-sample coupling is increased. At larger Tq an additional feature appears in the 
spectrum at about ±4 mV, which is due to transitions from the state \ip 2 ) to {ipa). 
These transitions are only possible because the probabilities pi are driven far from 
thermal equilibrium and p 2 has already a significant weight at \V\ =4 mV. Similar 
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bias V (mV) bias V (mV) 


Figure 12. Tunneling spectra of a spin S = 1 system at different coupling 
strengths Tq between tip and sample. The simulation parameters are identical 
to the ones used in figure [8] g = 2, D = —5 meV, E = 1 meV, f/ = 0, 
Jps = —0.1, T = 1 K. (a) Normalized zero field spectra for Tq ranging from 
1 X 10~^ (dashed red line) to 5 x 10~^ (full red line). At increased Tq a feature 
at eV = it I €3 — € 2 1 ~ ib4 meV appears and the apparent strength of the peaks at 
eV = dies increase, (b) State occupancy in steady-state for the smallest (dashed 
lines) and highest (full lines) coupling strengths, (c) and (d) Same as (a) and (b) 
but at an applied field of Bz = 10 T. The dashed-dotted line in (c) is the relative 
difference between the spectra at highest and lowest coupling strength. The inset 
in (c) shows the evaluation of the state energy ei with field. 

current induced pumping to higher excitation states has been observed for example 
for Mn-dimers adsorbed on CU 2 N [13], for small Fe clusters containing only a few 
atoms on Cu(lll) |2], or for Fe-OEP-Cl (Fe-octaethylporphyrin-chloride) molecules 
adsorbed on Pb(lll) [115j . Interestingly, the latter experiment was performed on a 
superconducting surface and with a superconducting Pb-tip which made it necessary 
to account for the gap around the Fermi energy in the quasi-particle density of states 
of tip and sample. 

If we now add a magnetic held of Bj, = 10 T to our simulation , the situation 
changes: The energy differences £12 and £23 between the hrst and second, and 
second and third eigenstate, respectively, are almost identical (see hgure inset) 
masking the additional step-like feature when pumping into energetically higher states. 
Nevertheless, also here a clear change of the spectrum occurs at high tip-sample 
couplings. We observe an increased relative conductance between £2 and £3 and a 
stronger third order peak at eV = ±£ 3 . This feature is mainly due to the third order 
scattering process (232) now possible, as apparently visible when taking the relative 
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(a) (b) (c) 




Figure 13. Simulated tunneling spectra of a Mn atom on Cu 2 N with spin S = 5/2 
at different coupling strengths Tq between tip and sample and at a magnetic field 
of Bz = 7 T. Simulation parameters are g = 1.9, D = —40 //eV, E = 0^ U = 1.3, 
Jps = —0.02, T = 1 K, and a tip polarization of r/* = 0.3. (a) Normalized spectra, 
and (b) averaged magnetic moment (rriz) in units of h for Tq ranging from 5x10“® 
to 2 X 10“^^, corresponding to a tunneling resistance of Rt ~ 270 to 6.8 Mfl 
or a stabilization current of 7 35 pA - 1.5 nA at a bias of eV = +10 meV. 

The dashed lines in (a) and (b) are calculated with 1/ = 0 at Tq = 2 x 10“^. (c) 
Occupancy of the six eigenstates at steady-state conditions for Tq = 2 x 10“^. 


difference between the spectra at strongest and weakest coupling strength (dashed 
dotted line in figure [T^t). 

As long as the tip is not spin-polarized, the current induced pumping into higher 
states cannot lead to an inversion of the state occupancy. The state’s lifetime for the 
eigenstate is inversely proportional to the sum of all scattering processes which 
leaves the state: 


rT = E + ry‘(er) + ry’^(eV)) . (42) 

3¥=i 

While for spin-averaged electrodes the scattering matrix elements do not change when 
inverting the initial and final state, the energetically higher states must have always 
a shorter lifetime. 

This behavior changes drastically when a spin-polarized tip is used. We have seen 
that, in particular for spin systems with a potential scattering U = S/2, the tunneling 
conductance and thereby the scattering rate depends strongly on the bias polarity 
(Figured]). Experimentally, this was first detected for Mn adatoms adsorbed on CU 2 N 
m and successfully discussed theoretically in a second order perturbative scattering 
model [SO] [79]. Figure [13] shows the simulated spectra of a Mn spin system at high 
field which - due to the small magnetic anisotropy - leads to an almost equidistant 
energy difference between the five eigenstates. These eigenstates are well described by 
pure states with the magnetic quantum numbers ruz = —5/2, —3/2,..., -1-5/2. 

At low tip sample couplings the dl/dV spectrum is similar to the ones simulated 
without rate-equations (Figure [9ji and e), but at higher coupling we observe a drastic 
reduction of the differential conductance at negative bias. This decrease is concomitant 
with the reduction of the average magnetization {rriz) = tr (^Sz J2iPi IV'i) (V’i/ ■ For 
the highest Tq the average magnetization becomes even negative at negative bias 
showing clearly the inversion of the state populations. Interestingly, the strong bias 
asymmetry in the dl/dV spectrum is only due to the potential scattering. Simulating 
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the system with U = 0 results in a much less asymmetric spectrum, even though the 
bias dependence of the state populations and average magnetization are unaltered. 

Finally, we have to remark that experimentally an effective interaction of the 
substrate conduction electrons with the Mn spin of Gs = 2.7 fiS was found m- 
This is much higher than the spin-substrate scattering determined solely by the 
spectroscopically estimated Jpo ~ —0.02 which results to Gs = 2e^jh x {JpsYS{S -\- 
1) « 0.3 /rS. This means that in this particular system roughly 90% of the spin 
relaxations with the substrate electrons must originate from scattering processes 
which do not directly leave their fingerprint in the observable differential conductance 
spectrum. In this context, ab-inito density functional calculations have found an 
about 3.1 times higher coupling to the substrate via the idxz and orbitals 

than experimentally observed |116j . while the approach discussed here neglect any 
orbital symmetry. 

5.1. Current induced correlations 

Comparing the dynamical model outlined above with experimental data that require 
the evaluation up to third order scattering is a very intriguing test of the capability, 
as well as the limitations, of this approach. Thus, we will now analyze spectroscopic 
data that have been measured for Fe adatoms on CU 2 N using spin-averaged and spin- 
polarized tips [HI]. Figure [HI shows the experimental data obtained for two different 
magnetic field directions; along the main anisotropy axis (easy axis, B^) which leads 
to a strong polarization of the eigenstates along the magnetic field (Figure [Tih l and, 
perpendicular to it, along the hard axis {Bx) which produces only a small polarization 
of the lowest eigenstates (Figure [T4|i). 

In this experiment the spin polarization was deliberately changed by vertical 
manipulation (’’picking-up“) [1171 [T^ of a Mn atom onto the apex of the tip, changing 
its polarization from rf 0 to rf ~ 0.4 and enabling the study of the identical atom 
with and without a polarized tip. The experimental data for the spin-average and 
for the spin-polarized measurements along the hard axis can be well simulated within 
one set of parameters (figure (Hb) . Surprisingly, the spectra are almost identical to 
simulations using only the zero-current approximation. Thus, the moderate coupling 
between tip and sample does not disturb significantly the state populations. This is 
also evident by plotting the average magnetic moment along the applied field (lower 
panels in figure [T4b and c). The small magnetic field (compared to the anisotropy 
energy) of only Bx = 3 T along the magnetically hard-axis cannot polarize the Fe 
atom significantly (Figure [T3]i). This means, that the spin imbalance produced by a 
spin-polarized current has only little influence on the dl/dV spectrum. However, note 
that the small shift towards higher energies of the step at \V\ ~ 4 mV, which has its 
origin in transitions between the states \tjji) and I'i/'a), might be due to an interaction 
between the spin system on the surface and the spin polarized electrode of the tip 

[nmiH]. 

The situation changes quite drastically when the magnetic field is applied along 
the easy axis and the spectrum is taken with the spin-polarized tip (Figure [Tib i. A 
strong peak at V Ri — 5 mV appears concomitant with the apparent disappearance 
of the step at precisely the same energy, that was clearly visible in the spectrum 
measured with the spin-averaging tip. Fortunately, the model describes this spectrum 
quite well, too and thus allows us to discuss the physical origin of these differences. 

Comparing the spectrum with the one calculated in the zero-current- 
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Figure 14. Correlations induced by tunneling substrate electrons, (a) Schematic 
state diagram for an Fe atom on Cu 2 N at Bz = 3 T. The blue arrows illustrate 
the third order transition (131). (b) and (c) Top panels: Experimental data from 
reference m of a Fe atom on CU 2 N measured with a spin averaging tip ( 77 * = 0) 
and a spin polarized tip (rf fs 0.4) at applied field in z (a) and x (b) direction 
(colored circles). The best fits (full lines) results in g = 2.11, D = —1.6 meV, 
E = 0.31 meV, U = 0.35, Jps = -0.1, Tgft = 1 K, and T§ = 3.2 x lO"'*. 
Spectra obtained and simulated with spin-polarized tips are shifted vertically by 
0.5 X 10~^{2e^/h). The dashed (dotted dashed) line in (b) is calculated with 
U = 0 (in zero-current approximation). Bottom panels: Corresponding average 
magnetic moment in field direction in units of h (full lines: spin-polarized tip, 
dashed lines: spin-average tip), (d) Schematic state diagram for an Fe atom at 
Bx = 3 T. (e) and (f) Schematic illustration (e) and Feyman-like diagram (f) of 
the third order scattering process (131) which leads to the peak at V Ri —5 mV 
in panel (a). 


approximation reveals that the disappearance of the step is indeed due to the increased 
coupling between tip and sample. The origin lies in the finite U, quite similar to 
the differential conductance reduction observed at negative bias for the Mn system 
(Figure [T31). Indeed, we see that the spectrum calculated with the same set of 
parameters except (7 = 0 leads to a spectrum in which the step-like increase of the 
dl/dV a,t V < —4 mV is maintained. Furthermore, we notice that the calculated 
magnetization (m^) decreased abruptly below this tunneling bias voltage, much faster 
than for the spin-averaging tip, due to the increased population of the states \ip 2 ) and 
1 ^ 3 ) (Figure [T3 }d lower panel). 

While the effects discussed above give a plausible explanation for the overall 
shape of the spectrum, there is still a noticeable discrepancy between experiment and 
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simulation. In particular, the experimental dl/dV data show a significantly stronger 
peak at ~ —5 mV than one would expect based on the simulation, resembling 
of the effects of strong correlations discussed in section U (which are not covered 
by the model). Thus, it is worth to dissect the characteristic peak which has its 
origin in the third order process (131) in which the spin system is scattered from 
|'0i) ll/’s) IV'i)- With the spin-polarized tip this transition has a significantly 
higher probability to occur when the electrons tunnel from sample to tip. The process 
is illustrated in figures [T4ji and[T4f. It starts with an j),) electron in the sample which 
scatters on the spin system changing its grounstate from |'0i) to I'i/'a) and then tunnels, 
with increased probability, as |t) electron into the majority states of the tip. The spin- 
down hole in the sample will lead to a slightly higher spin-up electron density close 
to the Fe atom. Next, a |t) electron from the sample fills the jj,) hole together with 
changing the localized spin-state from j^/^a) back to I'i/l’i) (Figure [Util. During this 
process the substrate electrons are correlated with the spin state. While the magnetic 
moment for \il)i) is approximately ruz ~ +2 and for j^/^a) Ri -|-1, respectively, at 
low enough temperature and at sufficiently negative bias and tunneling current, a 
correlated out-of-equilibrium Kondo-state develops which has the total wavefunction 

d/‘°‘“' = ^(|i^)|+2)-|f)|+l)). (43) 

Even though the experimental data only give hints that this correlated state has 
formed, a similar non-equilibrium Kondo formation has been found in transport 
measurements on carbon nanotubes |120j . The physical properties of such a non¬ 
equilibrium correlated state are very intriguing [12111122] . Entering this state, the 
substrate electrons partly screen the magnetic moment of the Ee spin reducing it to 
5—1/2. Such an underscreened Kondo state should show a very particular temperature 
and energy dependence [m nn mil ms]. Interestingly, in a system like the one 
discussed here, the transition between weak and strong coupling is not only governed 
by a characteristic temperature Tk |127j . but the formation of this correlated state 
is deliberately tuned by drastically increasing the probability of its creation when 
applying a spin-polarized current of sufficient energy by the probing tip. This enables 
us to envision for example pump-probe experiments m in which the time-evolution 
of the formation and the decay of this correlated state is measured in detail. 

6. Summary 

In this manuscript I have shown that applying perturbation theory to quantum spin 
systems enables one to describe experimentally measured differential conductance 
spectra with very high accuracy. This enables one to obtain a profound understanding 
of the physical processes on play and to separate single- as well as many-electron 
effects. 

The versatility of low-temperature scanning tunneling measurements on single 
and complexly coupled spin systems led me to believe that we should expect a 
multitude of exciting new experiments for the future, which will further deepening our 
fundamental knowledge on quantum systems in general and, in particular, quantum 
magnetism. Perturbative models, like the one outlined here, might be of help in such 
systems. Eor that, the supplemental material to this manuscript include an easy usable 
software package that allows not only to simulate the differential conductance spectra 
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of arbitrary complex spin systems, but additionally allows one to fit experimental data 
to the model. 

In this manuscript we have restricted ourselves to the adiabatic limit and 
neglected any time dependence in the parameters. However, it is straightforward 
to expand this model to capture also time dependent pump-probe measurements. In 
such a framework, coherent state superpositions could be accounted for by using, 
for example, a Bloch-Redfield approach in which the spin system is coupled to an 
open quantum-system |1281111411129] and in which interactions up to third order are 
included to account not only for the decay but also for the creation of coherences. 

Additionally, the study of the coupling between the spin-system and other (quasi)- 
particles, like phonons, photons, or magnons would be very interesting. Furthermore, 
while the perturbative approach fails when the system enters the strong-coupling 
Kondo-regime, it would be very intriguing to combine the simple model outlined here, 
with exact quantum impurity models. 
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Appendix A. Matrix elements for arbitrary spin polarization 


The two arbitrary spin density matrices p* and p® are a full description of the ensemble 
states in the tip and sample electron bath close to the Fermi energy. The eigenvectors 
\ip\) and |(/3|) of and p® are representative eigenstates of these incoherent ensembles 
and enables one to calculate the interaction transition intensities between them as: 




Aj, {‘Pf'l ^x,y,z,I ) 


(A.1) 


with A* and \\ as the eigenvalues of the corresponding eigenvectors. Note, that due 
to the in general incoherent spin ensembles in tip and sample equation I A. 1 1 has to be 
evaluated independently for all combinations oi i' = 1,2 and /' = 1,2. 

For the localized spin system we can write similarly the transition intensities 
between the eigenstates: 

= (A.2) 

so that the absolute square of the matrix element in equation [13] calculates to: 


M'mti 


=i:|Ke 








(A.3) 


The interference term between A4i and A42 (equation |23|) can be evaluated in an 
analogous way leading to: 

{x,y,z,I} 
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and for the reversed order processes: 

E eEs'.xe/'S^/x|?,^,5L-(A.5) 

j,k,l= i' m'f 

{x,y,zj} 

The tilde (^) and bar account for processes which starts and ends in the 

sample, and for processes in which tunneling electrons traverse the junction in the 
opposite direction, i. e. which starts in the sample and ends in the tip, respectively. 
They can be calculated in an analogous way with equation lA.ll 
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